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ABSTRACT 

We present a wide field census of resolved stellar populations in the northern half of M81, conducted 
with Suprime-Cam on the 8-m Subaru telescope and covering an area 0.3 square degrees. The 
resulting color-magnitude diagram reaches over one magnitude below the red giant branch (RGB) tip, 
allowing a detailed comparison between the young and old stellar spatial distributions. The surface 
density of stars with ages < 100 Myr is correlated with that of neutral hydrogen in a manner similar 
to the disk-averaged Kennicutt-Schmidt relation. We trace this correlation down to gas densities of 
~ 2 X 10^" cm~^, lower than typically probed with Ha flux. Both diffuse light and resolved RGB star 
counts show compelling evidence for a faint, extended structural component beyond the bright optical 
disk, with a much flatter surface brightness profile. The star counts allow us to probe this component 
to significantly fainter levels than is possible with the diffuse light alone. From the colors of its RGB 
stars, we estimate this component has a peak global metallicity [M/H] ^ —1.1 ± 0.3 at deprojected 
radii 32 - 44 kpc assuming an age of 10 Gyr and distance of 3.6 Mpc. The spatial distribution of 
its RGB stars follows a power-law surface density profile, I{r) oc r^'^ , with 7^2. If this component 
were separate from the bulge and from the bright optical disk, then it would contain ^ 10 — 15% 
of MSl's total V-band luminosity. Wc discuss the possibility that this is M81's halo or thick disk, 
and in particular highlight its similarities and differences with these components in the Milky Way. 
Other possibilities for its nature, such as a perturbed disk or the faint extension of the bulge, cannot 
be completely ruled out, though our data disfavor the latter. These observations add to the growing 
body of evidence for faint, complex extended structures beyond the bright disks of spiral galaxies. 

Subject headings: galaxies: individual (MSI) - galaxies: stellar content - galaxies: evolution - galaxies: 
structure 



1. INTRODUCTION 

The stellar outskirts of spiral galaxies hold vital clues 
to their formation and evolution. Under the framework 
of the ACDM cosmology, the halos of massive spiral 
galaxies arise larg ely from the merging and accretion of 
smaller sub-halos (|Bullock fc Johnston! I2005f ). The stel- 
lar populations of spiral galaxy halos are, thus, directly 
linked to the prope rties of the su b-halos and the accre- 
tion history (e.g., 'Bul lock fc Joh nston 2005; Font et all 
l2006f ). Furthermore, the merging and accretion that is 
especially common in ACDM at high redshift could be 
important in determining the properties of thick disks, 
which can dominate the total disk light at large radii 
and scale heights and which mounting evidence sug- 

*BASED ON DATA COLLECTED AT SUBARU TELESCOPE, 
WHICH IS OPERATED BY THE NATIONAL ASTRONOML 
GAL OBSERVATORY OF JAPAN. 
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Observing the outskirts of spiral galaxies is challeng- 
ing because these regions are very faint, typically at 
least several magnitudes below the sky level, or fxy ^ 
25 mag arcsec"^. Reliable surface brightness measure- 
ments at these faint levels require accurate flat-fielding, 
sky-determination, PSF modeling, and bright-star mask- 
ing. Diffuse light studies of halos and thick disks are 
usually limited to edge-on systems in which the thin disk 
light can be minimized by looking away from the mid- 
plane. For example, by stacking the images of over 1000 
edge-on galaxies in the Sloan Digital Sky Survey (SDSS), 
IZibetti et all (|2004f ) reached an r-band surface brightness 
/Zr ~ 31 mag arcsec"^. They found a halo of excess emis- 
sion around the stacked disk with an axis ratio of c/a 
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0.6 and projected intensity di stribution I(r) oc r , with 
7 = 2. In a related study, iZibetti fc FergusonI ()2004f ) 
examined the extraplanar hght distribution around a 
nearly edge-on disk galaxy in the Hubble Ultra Deep 
Field (HUDF) and found evidence for a halo with similar 
structural properties as the SDSS stacked halo. These 
structural properties ar e very similar t o those of the 
Milky Way (MW) halo ()Bell et al.l [20081 and references 
therein), suggesting that such halos are common around 
disk galaxies. However, the colors of the stacked SDSS 
and HUDF halos were much redder than expected for 
typical halo populations, suggesting an unusual stellar 
initial mass function (IMF), very hig h metallicity, or sig- 
nificant scattered light from the disk (jZibetti et al.ll20M : 
iZackrisson eTaI|[200l [de Jondl2008f) . 

Studies of the diffuse light around galaxies have also 
revealed evidence that thick disks are common (e.g., 
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tained optical and near-infrared images of many bulgeless 
edge-on disk galaxies reaching /i^ ~ 27 mag arcsec"^. 
From the color maps, vertical color gradients, and faint 
isophote shapes, they argued for the presence of red 
stellar envelopes around all their galaxies with surface 
brightnesses, spatial distributions, mean ages, and metal- 
licities similar to the MW thick disk. The unifor- 
mity in these properties is surprising if they were as- 
sembled from satellites with different masses, gas frac- 
tions, histories, and orbital properties, as expected 
in some hierarch i cal structure form ation models (e.g. 
lAbadi et al.]|2003l : iBrook et al.ll2004[) . However, the ages 
and metallicities were very uncertain because of the un- 
known star formation histories, the shallow depth of 
the near-infrared images, and the inherent difhculties 
of sky-subtracti ng and fiat-fielding in the infra red. Fur- 
ther analysis bv lYoachim fc DalcantonI ()2008bl ial) has in- 
deed revealed a variety of thick disk kinematics and 
compositions amongst the l ow-mass galaxies in the 
[Palcanton & Bernstein ()2002n sample. The properties 
of the higher-mass galaxies' thick disks were more dif- 
ficult to constrain because their thin disk light was too 
dominant in the regions with adequate signal. 

The best current method to study the faint outer 
regions of galaxies is to resolve their stellar popula- 
tions. This allows one to directly probe fainter sur- 
face brightness levels and to place tighter constraints 
on age and metallicit y than is pos si ble w ith diffuse 
light. For example, iMouhcine et al.l ()2005f ) studied a 
sample of nearby spiral galaxies with the Wide Field 
and Planetary Camera 2 (WFPC2) on board the Hub- 
ble Space Telescope (HST). They found the red giant 
branch (RGB) metallicity in the outskirts of these galax- 
ies, as inferred from the colors of their RGB stars, to 
be positively correlated with galaxy luminosity, qualita- 
tively consistent with theoretical expectations for stel- 
lar spheroids in ACDM (iFont e t al. 200l) and mono- 
lithic collapse scenarios (jLarsonl |1974| ). Their conclu- 
sions were complicated, however, by the small field-of- 
view of HST and the possibility of multiple populations 
along the line-of-sight. Indeed, when selected kinemat- 
ically, halo stars in the three L ocal Group spirals do 
not show a significant correlation ('Chapman et al .ll2006l: 
iKalirai e~al...2006. : .Ferguson et al..,2007; .Mouhcine et all 



[2007HMcConnachie et al"][2006l) . 

Ground-based, wide-field mapping of resolved stellar 
populations over large portions of nearby galaxies is an 
essential complement to pencil-beam HST-based studies. 
This approach is necessary for understanding their global 
structure, which may exhibit population gradients and 
inhomogeneities on a variety of scales. For examp l e, the 
INT and MegaCam surveys of M31 (llbata et al.l 120011 : 
[Ferguson et al.l 120021 : llbata et all I2007D have revealed a 
giant stellar stream and multiple other substructures dis- 
tributed throughout M31's halo which are undetectable 
with diffuse light and most easily seen in the spatial dis- 
tribution of RGB stars. These substructures have sizes 
and metallicity variations on scales much larger than 
HST's field-of-view. 

To place the most statistically significant constraints 
on theoretical galaxy formation models, we must increase 
the number of systems whose individual old and inter- 
mediate age stars have been resolved over large por- 
tions of their outer disks and halos. This means go- 
ing beyond the Local Group, a task which has only re- 
cently become possible thanks to the advent of wid e-field 
imagers on large tele scopes (e.g. Bland-Hawthor n et al.l 
120051: IDavid"gil2008allH ). At a di stance of 3.6 Mpc, cor - 
respondingto (m-M)o = 27.78 (|Freedman et al.lll994f) . 
M81 is the nearest massive spiral like the MW and M31, 
yet it is still close enough to resolve its low-mass giant 
stars. Wi th a Hubble morphological class of Sab and T- 
type of 2 (!de Vaucouleurs et al.lll99lT ). a total dynamical 
mass inside 20 kpc of ~ 10^^ Mp). and peak circular ro- 
tatio n velocity 250 km s""'^ (jCottesman fc Weliachewl 
|1975( ). M81 is similar in many respects to the MW and 
M31. It is therefore an essential system for testing our 
understanding of the formation of massive spirals. 

M81 is also the largest member of the nearest inter- 
acting galaxy group. Its two brightest neighbors, M82 
and NGC 3077, lie within a projected distance of 60 
kpc. The disturbed nature of the system is most easily 
visible when observing the 21 cm emission line of neu- 
tral hydrogen. Maps of the HI content throughout the 
system show extended tidal streams and deb r is between 
the g alaxies (jCottesman fc Weliachewl 
Il994f ). The computer simulation of 
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plains many of these tidal features as being the re- 
sult of close encounters between M81 and each of its 
neighbors ^ 200 — 300 Myr ago. The stellar contents 
of the most prominent debris, Arp's Loop and Holm- 
berg IX (HoIX ) hav e been studied with WFPC2 by 
iMakarova et al.l (l2002|) and with the Adva n ced Camera 
for Su rvevs /ACS) bvlde Mello et all (pOOl . Isn^bTeTan 
(|2008l ). and lWeisz et al.l (|2008l ). These authors concluded 
these objects may be tidal dwarf galaxies, new stellar sys- 
tems that formed in gas stripped from interacting galax- 
ies. 

A few papers have examined the large -scale, young 
resolv ed stellar content around M81 (e.g. 'D urrell et al.l 
120041 : iDavidge 2008b, 2009). They identified several 
groupings of 0-B stars and red supergiants which formed 
within the last ^ 100 Myr and which are located in M81's 
spiral arms and tidal tails. In this paper, we present 
wide-field imaging of M81 obtained with the Suprime- 
Cam instrument on the 8-m Subaru telescope. Our data 
reach over one magnitude below the RGB tip. This lim- 
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Fig. 1. — One square degree image showing the region observed with Suprime-Cam (large rectangle) relative to M81 and M82. North is 
up and east is left. The legend at upper left gives the scale of 10 arcmin (sj 10 kpc). Also shown are the locations of several other fields 
relevant to this paper which were imaged with Advanced Camer a fo r Surveys or Wide Field and Planetary Ca mera 2. These fields ar e 
lab eled 1-6 and were s tudied by (1) Tikhonov et al. (2003) and iMouhcine ct al. (2005), (21 ISabbi eFaP 1(2003 ) and IWeisz STan l|2008l 1. 
f Sl ide Mello et al.l ((20081 '). (4) and (5) de Jong ct al. (2008), and (6)[Williams ct al. (2008). 



iting depth is fainter than previous ground-based data 
in M81, and enables a detailed, yet global view of both 
young and old stellar structure around this galaxy. We 
present evidence for a faint extended structural compo- 
nent beyond the bright optical disk whose overall struc- 
ture and metallicity we can constrain because we resolve 
its stellar constituents. In [JH we outline the observa- 
tions and data reduction and, in we present the color- 
magnitude diagrams of point sources and non-M81 con- 
taminants. We examine the 2-dimensional stellar spatial 
distribution in 21 ^-i^d the radial surface brightness and 
stellar density profiles in We discuss possible inter- 
pretations in SJSland summarize the results in §7) 

Throughout this paper, we adopt for MBl's disk an 
inclination of 58 deg, and position angle (measured N 
through E) of 157 deg. At a distance of 3.6 Mpc, 1 arcmin 
« 1 kpc in projection. The radius of the galaxy is 
13.8 arcmin (|de Vaucouleurs et al.lll99j) . 

2. OBSERVATIONS AND DATA REDUCTION 

The obser vations were obtained w ith the Suprime-Cam 
instrument (jMivazaki et al.l I2002D on the 8-m Subaru 
telescope on the nights of January 7 - 8, 2005 (S04B, 
PI=N. Arimoto). This instrument consists of 10 CCDs 
of 2048x4096 pixels arranged in a 2x5 pattern, with a 
pixel scale of 0.2 arcsec and a total field of view of ap- 
proximately 34x27 arcmin (including long edge inter-chip 
gaps of 16 - 17 arcsec and short edge gaps of 5 - 6 arcsec). 

The original aim was to cover M81 using two field 



centers, one to the north and another to the south 
of its nucleus at {a, 5) = (9''55'"33^2, +69°03'55"). 
Weather conditions limited the amount of good-quality 
data obtained to the more northerly field centered at 
(9''55™30^0,+69°16'00") (J2000.0). The large rectan- 
gular box in Fig. [1] outlines the size and location of 
this field relative to M81, to M82, and to several other 
fields imaged with WFPC2 and ACS which are the 
most relevant to our current study of M81's stellar out- 
skirts. The HST fields a rc labeled 1 ~ 6 and were stud- 
ied by (1) iTikhonov eFal. (2 005) and IMouhcine et all 
imm. (2) ISabbi et all (120081) and IWeisz et all (l200aH 
(3) I de Mello et al.l (120081), (4 ) and (5) Ide Jong et all 
(20081). and (6;i lWiniams etaH jlOOS). 

We obtained a set of 10 images of M81 in the Johnson 
V filter with individual exposure times of 630s on each 
night, and 20 images in the Sloan i' filter with exposure 
times of 215s on the night of Jan. 7th. All observations 
were recorded under non-photometric conditions through 
patchy thin cirrus. Comparison of expected throughput 
with subsequent calibration based on SDSS photometry 
for this region indicated that the average attenuation for 
these images was < 10%. The set of 10 y-band images 
obtained on Jan. 7th in relatively poor seeing (~ 1.5") 
were later rejected from the analysis. The remaining V- 
band images were taken in an average seeing of '~ 0.7" 
while the i'-band images had an average seeing of ~ 1.1". 
To fill in the chip gaps and facilitate the removal of cos- 
mic rays and bad pixels, individual images were dithered 
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by ~ 25 arcsec, resulting in a mosaic with total field 
of view « 36x28 arcmin. Flat-field and interchip gain 
variations were removed with master flats obtained by 
combining 12 and 11 twilight sky flats in the V and i' 
filters, respectively. After flat-fielding, remaining large 
scale variations in dark sky level, measured directly from 
stacked dark sky images at several different positions ob- 
tained during this run, were less than 1% of sky. An 
j'-band fringe frame acquired from an earlier Suprime- 
Cam imaging run was used to help assess the degree of 
dark sky fringing present, but this was found to be neg- 
ligible in our data, so this extra image processing step 
was not required. 

After converting the raw data to multi-extension FITS 
format, all images and calibration frames were run 
through a variant of the data reduction pipeline de- 
veloped for the INT Wide Field Survey ^. The main 
steps we used in producing a detected object cata- 
logue, including background estimation, object detec- 
tion, parameter estimation, and mor phological classifca - 
tion, are descr i bed in more detafl inll7^(fl985l ll99l. 
Ilrwin fc Lewii (|2001h . and llrwin et al.l (|2004l ) . Here we 
present a brief overview of these steps. 

Prior to deep stacking, catalogues were generated for 
each individual processed science image to both refine the 
astrometric calibration and asses the data quality. For 
astrometric calibr ation, we found that a Zeni thal poly- 
nomial projection (jGreisen fc Calabrettall200§ ) provided 
a good prescription for the World Coordinate System 
(WCS) and included all the significant telescope radial 
field distortions. We used this in conjunction with a 6- 

http:/ /www. ast.cam.ac.uk/wfcsur/ 



parameter linear plate model per detector to define the 
required astrome tric tra nsformations. The 2MASS point 
source catalogue (|Cutri|[2003i) was used for the astromet- 
ric reference system. 

The individual image qualities were then assessed us- 
ing the average seeing and ellipticity of stellar images, 
as well as sky level and sky noise determined from the 
object cataloging stage. As a result of this assessment, 
the set of V-band exposures taken on Jan. 7th, in much 
poorer seeing conditions, were not included in the final 
V-band stack. Images were stacked at the detector level 
using the updated WCS information to accurately align 
them with respect to a reference image. The background 
level in the overlap area between each stacked image and 
the reference was adjusted additively to compensate for 
sky variations during the exposure sequence and the fi- 
nal stack included seeing weighting, confidence (i.e. vari- 
ance) map weighting, and clipping of cosmic rays. 

Next, we generated detector-level catalogues from the 
stacked images and updated the WCS astrometry in the 
FITS image extensions prior to mosaicing all detectors 
together.^ Residual astrometric errors over the whole 
stacked array were typically < 0.1", greatly simplifying 
this process. Slight offsets in underlying sky level be- 
tween the stacked detector images caused small (typi- 
cally '^0.1-0.2% of sky), but still visible, discontinuities 
in the final mosaics. These offsets were due to small color 
equation differences in the detectors and the relatively 
blue color of the twilight sky compared to dark sky and 
unresolved diffuse light from M81. We corrected these 

^ Mosaicing the ensemble to form a single large image minimizes 
the impact of lower total effective exposure times in the overlap 
regions between detectors, particularly toward their edges. 
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Fig. 3. — Color-magnitude diagram (left) and Hess diagram (right) of ~ 40,000 point sources in the Suprime-Cam mosaics of M81. The 
Hess diagram is a 2-dimensional histogram showing source density on a logarithmic intensity scale. Error bars on the right-hand side of 
each panel indicate the median photometric error at each magnitude. 



offsets iteratively by visual inspection of a 4x4 blocked 
version of the mosaics using a pre-assigned keyword in 
each relevant detector FITS extension designed for this 
purpose. 

Fig. [21 shows a color composite image constructed from 
the V and i' mosaics, where north is up and east is to the 
left. The image has a linea r intensity scaUng an d a color 
mapping similar to that of iLupton et al.l (|2004D with V 
for the blue channel, i' for red, and the average for green. 
The edges were cropped due to the slightly different field 
coverage in the final stacked mosaics. 

The background light variation over these mosaics was 
rather complex and varied on relatively short angular 
scales. Therefore, as an alternative to our default method 
of background estimation and removal, we also used an 
independent method prior to cataloging in the crowded 
central region of M81. This step used a filter based on 
an iterative, non-linear unsharp masking and clipping 
algorithm to progressively remove fine details to a user- 
chosen scale length. An exa mple of the benefits of thi s 
approach is shown in Fig. 1 of lMcConnachie et al.l (|2007f ) , 
where complex background variations due, in that case, 
to scattered light were effectively removed. Comparison 
with photometry generated in a more conventional way 
demonstrated that even for the brightest unsaturated 
stars this aggressive filtering made less than 1% differ- 
ence to their measured fluxes and had completely negli- 
gible effect on fainter stars. The significant advantages 
were that if this filtering step was run prior to mosaicing 
it obviated the need for detector-level adjustments and 
it greatly simplified subsequent object detection and flux 
estimation. 

Object fluxes were estimated using "soft-edged" aper- 
tures of radius w FWHM, which, for relatively un- 
crowded images, deliver 80-90% of the signal-to-noise of 
a perfect PSF fit in a fraction of the time (e.g., Irwin 
ll997HNavloilll998D . A significant portion of the detected 
objects in these images were background galaxies, for 
which aperture photometry was more useful than PSF- 
fitting. Therefore, a series of apertures ranging from 1/2 
to 4 times the FWHM were additionally used as the basis 



for the morphological classification in a curve-of-growth 
analysis. We also used these apertures to compute stellar 
aperture corrections for the primary fiux measure. The 
PSF varies in a rather subtle way over such wide-field 
images, and even though aperture photometry is not as 
sensitive to this variation as PSF-fitting, we still used the 
different aperture fluxes to check and correct for it. 

Beyond M81's bright optical disk, the typical object 
separation increases from ~ 5" at a deprojected radius 
Rdp = 14 arcmin to ^ 10" at 40 arcmin. To mitigate 
the effects of crowding, which are more important at 
Rdp < 14 arcmin, the pipeline employed two extra 
refinements to the flux measurement. First, the object 
detection information for neighboring regions was used to 
flag those pixels corrupted by external groups of detected 
objects. Second, object fluxes from a pixel image blend 
of overlapping sources were solved for simultaneously via 
least-squares fltting of the aperture fluxes, i.e. equivalent 
to a top-hat shaped PSF, using flxed input coordinates. 
With this approach, bad pixels and regions of extremely 
low confldence were also readily flagged and avoided ap- 
propriately. Direct comparison of the pipeline-generated 
photometry with full PSF-fltting over a broad range of 
images for other projects has shown that in all but the 
most densely crowded regions, like the centers of globular 
clusters, similar quality results are achieved. 

Due to the non-photometric observing conditions, we 
brought our flux measurements onto a Vega-like system 
(Johnson V and Cousins /) by matching stars in our cat- 
alog with those in the SDSS, after applying the Lupton 
transformation equations provided on the SDSS website 
. Only stars with photometric errors in the transformed 
system < 0.2 mag were used in the matching. The re- 
sulting photometric zero-point is accurate to ~ 0.05 mag, 
which is adequate for our study. 

The final catalog contains ~ 40,000 objects with stel- 
lar, or probable stellar, classifications in both passbands 
(hereafter referred to as point sources), meaning their im- 
age shape parameters lie within 2 — 3(t of the main stellar 

^ http:/ /www. sdss.org/dr6/algorithms/sdssUBVRITransform. html 
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locus. Relaxing this criterion in / makes little difference 
to our results, but relaxing it in V actually increases con- 
tamination from unresolved background galaxies since 
the seeing in V was significantly better. The morpholog- 
ical classification algorithm categorized ^ 23, 000 objects 
as background galaxies, which we discuss in more detail 
below. 

3. COLOR-MAGNITUDE DIAGRAM 

Fig. [3] shows the final point source catalog as a color- 
magnitude diagram (CMD) in the left panel and as a 
Hess diagram in the right panel. The Hess diagram is a 
2-dimensional histogram showing the source density on 
a logarithmic intensity scaling. There are several clear 
stellar sequences visible, which we discuss in more detail 
below. On the right-hand side of each panel are the me- 
dian photometric errors reported by the data reduction 
pipeline for all point sources in 0.5 mag- wide bins, where 
the color errors are the V- and /-band errors added in 
quadrature. The median J-band error is 0.1 at / ~ 24.7 
and w 0.14 at / = 25 while the median color error is 0.1 
at / w 24.3 and « 0.16 at / = 25.0. These errors reflect 
Poisson noise in the sources and the sky level. They do 
not account for systematic color and magnitude shifts or 
for correlated errors in both bands, which can occur in 
heavily crowded regions, but by restricting most of the 
point source analysis to M81's outskirts, we do not ex- 
pect this to be a serious issue. We also note that point 
sources are saturated at V < 19 and / < 19. 

The extinction maps of lSchlegel et al.l (| 1998( 1 indicate a 
variable extinction across our field, with ^ 15% of point 
sources having E{B — V) values at least 0.01 mag in 
excess of the median (0.08) and a total range of 0.05 
- 0.12. Henceforth, we apply extinc tion co rrections on 
a st ar-by-star basi s using thc'SchlcgefeEaD ^98) maps 
and lCardelli et ahl (|1989( ) extinction law, for which Ry = 
3.1 and Ai/Ay = 0.479. ISchleeel et all (|1998f) removed 
extragalactic objects from their maps down to a certain 
flux limit and replaced them with median values from 
the surrounding sky. By examining their mask maps, one 
can see that the inner ^ 16 arcmin in projected radius 
around M81 (including HoIX and part of Arp's Loop) 
was indeed processed in this manner. Thus, the star- 
by-star correction does not include extinction internal 
to M81, but this should not be a serious problem since 
wc focus on the RGB stars, which should be the least 
affected by internal extinction (Zaritsky 1999), and our 
CMD selection boxes (described below) are larger than 
the expected amount of internal extinction. For instance, 
the spectro-photom etric study of M81's diffuse light by 
iKong et all (pOOOl ) found an average E{B — V) ^ 0.2 
throughout M81's bright optical disk, suggesting there 
may be ^ 0.12 mag of internal reddening within this 
region. 

In Fig. 01 we show the p oint source H e ss dia gram 
overlaid with isochrones from iMarigo et al.l (|2008[ ) with 
the circumstellar dust option turned off. The young 
isochrones at (^ — /) ~ have ages of 10.0, 17.8, 
31.6, 56.2, 100, and 178 Myr and a metallicity [M/H] = 
\o%{Z I Z(£) — —0.4. We use this metallicity as a fidu- 
cial value because the young stars sampled by our field 
exist in a variety o f environments with JM/H] ranging 
from ry -0.7 to fZaritsky et al."1994'; Ide Mello et all 
[2008t iSabbi et al.|[2008; .Davidge.2008bl . The three old 
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Fig. 4. — Hess diagr am of ~ 40, OOP poi nt sources with theo- 
retical isochrones from IMarigo et al.l II2008I) overlaid. The young 
isochrones at (V — /) ~ have ages of 10.0, 17.8, 31.6, 56.2, 
100, and 178 Myr and a metalUcity [M/H] = -0.4. The three 
old isochrones at (V — /) ~ 1 — 3 have a common age of 10 Gyr 
and [M/H] = -1.3, -0.7, and -0.4. The boxes are used to select 
M81 stars in different evolutionary stages: main sequence and blue 
helium burning (MS+BHeB), red supergiant (RSG), asymptotic 
giant branch (AGB), and red giant branch (RGB). The 50% com- 
pleteness level is estimated to lie at /q ~ 25.0 for MS-|-BHeB stars 
and Iq ~ 24.4 for RGB stars (see [[3]for details). 



isochrones at (y — /) 1 — 3 have a common age of 10 
Gyr and [M/H] = -1.3, -0.7, -0.4. The discontinuities 
in the asymp t otic g iant branch (AGB) are explained in 
IMarigo et al.l (|2008[ ) and are caused by changes in the 
opacity tables at the transition to the thermally pulsing 
phase. 

In what follows, we will focus on several particular 
CMD regions (outlined in Fig. which isolate stars in 
different evolutionary stages at M81's distance. The blue 
lines in Fig.|4|mark the region occupied by main sequence 
stars < 20 Myr old and by blue helium burning stars at 
the hottest extension of the blue loop phase (MS-fBHeB) 
with ages < 100 Myr. Stars within the cyan polygon are 
red supergiants (RSGs) with ages in the range ~ 20 — 200 
Myr. The tip of the RGB lies at / ~ 24, so stars within 
the red lines are mainly RGB stars with ages ^1 — 10 
Gyr. Note that there could be some contamination of the 
RGB box by young, red helium burning stars with masses 
of - 3 - 4 Mq, particularly if they have [M/H] > -0.4. 
The magenta lines enclose AGB stars above the RGB tip, 
which tend to have some what younger ages (^ . 5 — 8 
Gyr) than the RGB stars (jMartmez-Delgado et aI1ll999l : 
iGallart et al. 2005). 

We have also divided the RGB box into "metal-poor" 
([M/H] < -0.7) and "metal-rich" ([M/H] > -0.7) sub- 
regions. There could be some overlap in the metallicities 
probed by these subregions, but they are broader in color 
than the photometric errors, and so are useful in identi- 
fying any population gradients. Some of the most metal- 
poor ([M/H] < -1.3) and metal-rich ([M/H] > -0.4) 
RGB stars may fall outside the total RGB box, but ex- 
tending it further to the blue or red would increase con- 
tamination from MW foreground stars and unresolved 
background galaxies, and increase uncertainties due to 
incompleteness. Comparing the number of RGB stars 
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Fig. 5. — Hess diagram o f foreground stars b ased on the Besangon 
model of the Milky Way (I Robin et all 120031 ') for a field with the 
same area and line of sight as the M81 mosaics. The stars were 
scattered according to the photometric errors of the real data, but 
no completeness corrections were applied. The number of predicted 
foreground stars is ~ 2% of the number of point sources over the 
observed magnitude range. The selection boxes (from Fig.lJJ avoid 
the most heavily contaminated regions. 



with [M/H] > -0.4 in the iWilliams et all (pOOl CMD 
(field 6 in Fig. [T] at R^p = 14 arcmin) to the number in 
our Subaru CMD reveals that most of these stars are too 
faint to be detected in our data. However, our detailed 
RGB analysis focuses on exterior regions where signifi- 
cantly fewer metal-rich RGB stars are expected to exist, 
a point substantiated in |j5l 

As can be seen in Fig. E] - [3 these selection boxes 
sample regions of the CMD that are likely to maximize 
the number of M81 stars relative to the number of fore- 
ground MW stars and background galaxies. Fig.Oshows 
the foreground star Hess diagrani predicted by the Be- 
sangon model ([Robin et al.l l2003f) for a field with the 
same area and line of sight as ours. The stellar colors and 
magnitudes have been scattered using a simple exponen- 
tial function to mimic the increase of photometric error 
with magnitude seen in the real data. The Besangon 
model predicts that the number of foreground stars is 
~ 2% of all point sources over the magnitude limits of 
the observed CMD. Since applying completeness correc- 
tions would only decrease this already small percentage, 
we do not use this model to correct the stellar number 
counts for MW stars. 

At bright magnitudes (/ < 23), morphological classifi- 
cation effectively removes background galaxies, but they 
can still be a significant source of contamination at faint 
magnitudes, especially when observing from the ground. 
To check the effectiveness of our classification algorithm, 
Fig. [5] shows the Hess diagram of the ~ 23,000 objects 
classified as extended. The CMD selection boxes are 
overlaid to facilitate comparison with the other CMDs 
presented here. The extended objects are concentrated 
in a broad diagonal band, the bulk of which lies at bluer 
colors than the RGB and AGB selection boxes. Impor- 
tantly, the extended objects are distributed differently 
from the point sources, with no obvious stellar sequences 
aside from some misclassified MS-|-BHeB stars located in 
spatially crowded regions where their PSFs overlap sig- 
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Fig. 6. — Hess diagram of ~ 23, 000 extended objects in the M81 
field. The selection boxes (from Fig. |4]l avoid the most heavily 
contaminated region. The surface density of extended objects in 
the total red giant branch box is less than half that of the point 
sources in Fig.|4]at all radii. 



nificantly. This fact suggests that the algorithm has cor- 
rectly classified the majority of detected M81 stars. The 
surface density of point sources within the total RGB box 
is over 2 times higher at all radii than the surface density 
of extended objects in the same box. Nevertheless, some 
compact background galaxies can be mistaken for point 
sources and these contaminants become increasingly im- 
portant at large galactocentric distances. 

The large spatial extent of M81 and various tidal fea- 
tures throughout Suprime-Cam's field-of-view make it 
difficult to estimate the contamination from background 
galaxies misclassified as point sources using just the M81 
mosaics themselves. However, we can estimate this level 
from images of the edge-on spiral galaxy, NGC 4244, 
taken during the same observing run. Fig. [7] shows the 
Hess diagram of a control field extracted from an area 
of 112 arcmin^ at vertical heights > 20 kpc above NGC 
4244's disk, assuming a distance of 4.4 Mpc (jSeth et al.l 
120051 ). This area corresponds to ~ 10% the total field-of- 
view around M81. Extinction corre ctions were applied to 
each object individually using the lSchlegel et all (|1998D 
maps, which indicate a median E{B — V) of 0.02. The 
control field has ~ 2, 000 sources distributed in the CMD 
in a similar manner as the extended objects of the M81 
field. The most densely populated region lies in the same 
diagonal band at0<(F-/)o<l. The CMD selection 
boxes avoid this region while still sampling as much of 
M81's stellar populations as possible. The surface densi- 
ties of objects within the boxes in the contaminant CMD 
are listed in Table [1] with Poisson errors. The contami- 
nant surface density in the total RGB box is less than the 
corresponding point source density at all radii, reaching 
a maximum fraction of ~ 50% at a deprojected radius 
Rdp 37 arcmin. 

We also estimated background contamination by 
downloading from the HST archive two M81 fields ob- 
served with ACS (SNAP 10523, PI: de Jong) at Rdp ~ 
20 — 30 arcmin. These two fields were observed in the 
F606W^ and F814W^ filters and are indicated as fields 4 
and 5 in Fig. [1] We performed the data reduction using 



8 



Barker et al. 



^ 22 




Fig. 7. — Hess diagram of 2, 000 contaminants extracted from 
an area 112 arcmin^ 10% the total field-of-view around M81) 
at vertical heights > 20 kpc from the plane of NGC 4244. These 
objects are mostly unresolved background galaxies. The selection 
boxes (from Fig. Q avoid the most heavily contaminated region. 
Table [l] lists the number counts of objects in the boxes, which 
are used to subtract contaminants from the radial stellar density 
profiles of M81. 



TABLE 1 
Surface densities of contaminants. 

CMD NGC4244 ACS 

region (arcmin~^) (arcmin~^) 



MS+BHeB 

RSG 

AGB 

Metal-poor RGB 
Metal-rich RGB 
Total RGB 



0.03 ± 0.02 
0.37 ±0.06 
0.52 ±0.07 
1.28 ±0.11 
0.64 ±0.08 
1.93 ±0.13 



0.3 ±0.2 
0.3 ±0.2 
0.4 ±0.2 
2.0 ±0.4 
1.0 ±0.3 
3.0 ±0.5 



the ACS module of the DOLPHOT package " following 
the basic steps outlined in the DOLPHOT manual. We 
defined objects as point sources if they were classified by 
DOLPHOT as 'good stars' with S/N > 5 and crowd- 
ing parameter < 0.5 in both filters, and if the overall 
\sharp\ < 0.1 and x < 3. Magnitudes in the ACS filter 
system were transfor med to the Johnson- C ousins system 
using the equations in lSirianni et al.l (|2005( ). We matched 
point sources in our Subaru catalog with point sources 
in the ACS catalogs by applying small, constant off- 
sets in right ascension and declination and then counted 
the number of unmatched point sources in our catalog, 
presumed to be unresolved background galaxies. Note 
that this method assumes the ACS point-source catalogs 
themselves have no contamination from unresolved back- 
ground galaxies. Reassuringly, the resulting background 
surface densities in each of the CMD selection boxes (see 
Table [T|) were consistent with those derived from the 
NGC 4244 data to within ~ 2a. Compared to the NGC 
4244 background, the ACS background is ^ 50% higher 
in the total RGB box, but the ratio of metal-poor/metal- 
rich RGB background is the same. The ACS background 

DOLPHOT is an adaptation of the photometry pack- 
age HSTphot (Dolphin 2000). It can be downloaded from 
[http: / /purcell. as. arizona.edu/dolphot/ 
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Fig. 8. — The dereddened luminosity function of point sources 
in two combined ACS fields (open histogram; fields 4 and 5 in 
Fig. [T} compared to that of the matched point sources in our cata- 
log (hatched histogram). The ratio of the two luminosity functions 
indicates that our catalog is > 50% complete for Ig < 24.4 at 
deprojected radii R^p = 20 — 30 arcmin. 



is also higher in the MS+BHeB box, but lower in the 
AGB and RSG boxes. 

We decided to use the NGC 4244 data to subtract the 
background from the radial density profiles presented be- 
low, since they sample an area 5 times greater than the 
ACS fields, are free from the uncertainties in the object 
matching, and rely on the same instrument, filters, data 
reduction procedure, and morphological classification as 
the M81 data presented here. This approach assumes a 
negligible contribution from stars and globular clusters 
in NGC 4244's halo 20 kpc above the plane and that the 
background galaxy and foreground star populations (at 
the magnitudes we are probing) are statistically equiv- 
alent to those around M81. At / ~ 23 — 25, the back- 
ground galaxies have a median z 1 (lllbert et al. 2006) . 
At this redshift, 1" « 8 kpc, and the entire Suprime- 
Cam field-of-view samples a projected volume ~ 15 Mpc 
wide, large enough to gather a statistically representa- 
tive sample of background galaxies regardless of position 
on the sky. Cosmic variance may be more of a con- 
cern here, since we are not able to use Suprime-Cam's 
entire field-of-view. However, cosmic variance is more 
likely to affect the ACS background, which covers much 
less area. With regards to foreground contamination, 
NGC 4244 is separated from M81 by « 40 deg and lies 
at a higher Galactic latitude, but both galaxies are sig- 
nificantly above the MW plane. Indeed, the Besangon 
model predicts only ~ 10 — 25% fewer foreground stars 
towards NGC 4244 than towards M81, within the CMD 
selection boxes. Our general conclusions are unchanged 
if we adopt the ACS background estimates, but to be 
conservative, we adopt the differences between the ACS 
and NGC 4244 estimates as systematic errors, which are 
added in quadrature with the Poisson noise in the latter. 

Another important issue is estimating the complete- 
ness rate of our Subaru catalog. The standard method 
is to run artificial star tests, in which thousands of stars 
with known magnitudes and colors are inserted into the 
original images, which are then re-processed. However, 
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TABLE 2 
Luminosity function 

PEAK magnitudes. 



(V - I)o LF peak 
range Iq magnitude 



-1.0 - 0.5 25.15 

0.5 - 1.0 24.85 

1.0 - 1.5 24.55 

1.5 - 2.0 24.45 

2.0 - 2.5 24.05 

2.5 - 3.0 23.35 



this method requires detailed knowledge of the PSF and 
how it varies across the field-of-view. Since we have 
used optimized aperture photometry that is insensitive 
to these variations, and given the complexities of the mo- 
saicing and data reduction, we have elected not to run 
artificial star tests. Instead, we have used two alternative 
methods to probe the completeness rate of our catalog. 

First, the completeness rate can be estimated by com- 
paring the number of point sources in the Subaru and 
ACS catalogs. Fig. [8] shows the dereddened luminosity 
functions of all point sources in the two ACS background 
fields (open histogram) and the matched point sources 
in the Subaru catalog (hatched histogram). This figure 
shows that in the range Rdp ~ 20 — 30 arcmin and for 
7o < 24.4, the Subaru catalog is > 50% complete relative 
to the ACS catalog. In addition, the number counts in 
the CMD selection boxes indicate that each box is > 50% 
complete relative to the ACS catalog, with a statistical 
uncertainty of 5% or more depending on the number of 
sources in each box. 

Second, we use the Subaru luminosity function (LF) 
peak as a proxy for the 50% completeness level. It is only 
an approximate indicator of this level because it depends 
on the true shape of the LF. Nevertheless, for our entire 
catalog, this peak occurs at Iq ~ 24.45, consistent with 
the above estimate. Table [5] lists the LF peak magnitude 
for several different color ranges. This table shows that 
all of the MS-fBHeB and RSG boxes and most of the 
AGB box lie brighter than the LF peak, while about half 
the metal-poor RGB box and most of the metal-rich RGB 
box lie fainter than the peak. We note that restricting 
our analysis to stars brighter than the LF peak makes no 
significant difference to our conclusions, but we use the 
entire CMD boxes, anyway, to decrease Poisson noise. 
Thus, one should be mindful of the different completeness 
levels for the metal-poor and metal-rich RGB stars when 
examining their relative distributions. 

4. TWO-DIMENSIONAL SPATIAL DISTRIBUTION 

Fig. [5] displays the tangent plane projection of all point 
sources with Iq < 25.15 and Vq < 26.85. The ellipses cor- 
respond to deprojected radii of Rdp = 10, 20, 30, 40, and 
50 arcmin (recall that 1 arcmin « 1 kpc). No correc- 
tion for contaminants has been made. The diffraction 
spikes of a few highly saturated stars appear as thin ver- 
tical white stripes. Visual inspection reveals significant 
substructure associated with spiral arms in M81's disk 
and the tidal features, HoIX at (^,7?) - (0.2, -0.02) and 
Arp's Loop at ~ (0.15,0.25). The hole in the central 
region is due to saturation and stellar crowding. 




a.3 0.2 0.1 0.0 -0.1 -0.2 -0.3 

Fig. 9. — Tangent plane projection of point sources with Iq < 
25.15 and Vq < 26.85. Ellipses denote deprojected radii of 10, 
20, 30, 40, and 50 arcmin (1 arcmin fsi 1 kpc). No correction for 
contaminants has been made. The diffraction spikes of a few highly 
saturated stars appear as thin vertical white stripes. 



Going clockwise from top left. Fig. [TOl shows the spa- 
tial maps of RGB, AGB, RSG, and MS+BHcB stars. 
No correction for contaminants has been made. From 
these maps, one can see that most of the overdensities 
observed in Fig. [H] are due to star formation in the last 
100 Myr. There are no stellar groupings in the RGB and 
AGB maps that can be unequivocally linked to group- 
ings in the MS+BHcB or RSG maps. Indeed, the RGB 
and AGB stellar distributions are smooth except near the 
chip edges and mosaic corners, where the completeness 
rate is lower due to different field centers and vignetting. 
These regions are therefore excluded from the radial stel- 
lar density profiles discussed in [JH 

In their ACS imag es of Arp's Loop at (0.16,0.23), 
Ide Mello et"all (|2008i f found that the RGB stars were 
broadly distributed like the youngest stars, except for 
a couple O-B associations, which may be too young to 
have formed RGB stars. Our stellar maps, which reveal 
the stellar distribution on a larger scale than the HST 
images, appear to contradict this finding because there 
is no clear broad overdensity of RGB stars associated 
with Arp's Loop. This fact is consistent with the idea 
that Arp's Loop was initially entirely gaseous tidal debris 
that began forming stars ^ 200 — 300 Myr ago, around 
the time of the last close encounters between M81, M82, 
and NGC 3077. It is possible that young, red He-burning 
stars were contaminating the Ide Mello et all jlOOS) RGB 
sample, but there is no obvious reason why their sample 
would suffer more from this contamination than ours. It 
is also possible that the completeness rate in our maps is 
not exactly uniform across this part of Arp's Loop. This 
is obviously the case for HoIX, where the crowding is so 
severe that it appears as a hole in the RGB star map, so 
we are unable to comment on the assertion of lSabbi et al.l 
(|2008[ ) that most of its RGB stars belong to M8L 

One o f the stellar concentrations identified bv lDavidgi 
(|2008b( ). TDO 3, lies within our field-of-view at (- 
0.13,0.25) and appears as an overdensity of MS+BHeB 
and RSG stars. TDO 3 is clearly associated with the 
northernmost sp iral arm, which may be tidally driven 
(|Yun et al.lll994[) . but neither of these objects necessar- 
ily are, or will become, distinct tidal dwarf galaxies, as 
may be the case for Arp's Loop and HoIX. The CMD of 



10 



Barker et al. 




0.3 0.2 0.1 0.0 -0.1 -0.2 -0.3 




0.3 0.2 0.1 0.0 -0.1 -0.2 -0.3 




Fig. 10. — Going clockwise from top left, the tangent plane projection of red giant branch stars (ages 1 — 10 Gyr), asymptotic giant 
branch stars (ages ~ 0.5 — 8 Gyr), red supergiant stars (ages ~ 20 — 200 Myr), and main squence and blue helium burning stars (ages 
< 100 Myr). Ellipses denote deprojected radii of 10, 20, 30, 40, and 50 arcmin (1 arcmin x 1 kpc). No correction for contaminants has 
been made. 
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Fig. 11. — Tangent plane projection of main sequence and blue 
helium burning stars (ages < 100 Myr). Ellipses denote depro- 
jected radii of 10, 20, 30, 40, and 50 arcmin (1 arcmin 1 kpc). 
No correction for contaminants has been m ade. Contours sho w the 
column density of HI in the VLA mosaic of lYun et al.l II1994I ). The 
contour levels increase by factors of two with the lowest (white) 
covering < 2 X lO^'' cm"^ (or 1.7 Mq pc"^). 

TDO 3 suggests an age spread of several tens of Myr with 
the young est stars ~ 20 M yr old, consistent with the es- 
timates of iDavidgd (|2008bD . We are also able to confirm 
his tentative detection of a stellar grouping east of Arp's 
Loop at (0.25,0.24), with a similar age as TDO 3. 

Fig.rrrishows the distribution of MS +BHeB stars over- 
laid on the HI column density map of lYun et~al., (,19941 . 



The contour levels increase by factors of two with the 
lowest (white) covering densities < 2 x 10^° cm~^ (or 
1.7 Mq pc~^). The young stars closely trace the HI for 
column densities > 8 x 10^° cm~^. Few young stars exist 
in regions with lower column density. This seems consis- 
tent with a simple star formation threshold gas surface 
density, but it could also be caused by a star formation 
rate in these low-density regions that is too small to pro- 
duce the high-mass stars we have resolved. 

To investigate the correlation between HI and 
MS+BHcB stars in more detail, we divided the field into 
square bins 1.7 arcmin (w 1.7 kpc) on a side and counted 
the number of these stars in each bin. For reference, 
the beam size of the radio observations was ~ 1 arcmin. 
Fig. [T2] shows the background corrected number counts 
as a function of mean HI column density within each 
region. Errors with downward pointing arrows indicate 
regions with only one young star. A maximum likelihood 
fit yields a best fit slope of 1.6 and 2a confidence inter- 
val of ±0.2. Increasing or decreasing the area of each 
bin by a factor of 4 increases or decreases the slope by 
~ 0.1 — 0.2, respectively. 

Our measured slope is roughly consistent with expec- 
tations from the canonical Kenni cutt-Schmidt (K S) re- 
lation, which has a slope of 1.4 ([Kennicu tt HToM ). but 
that relation was derived for disk-averaged quantities 
rather than individual regions. Our slope is within 
the range (1.2 - 3.5) derived for small regions or az- 
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Fig. 12. — Background-subtracted surface density of main se- 
quence and blue helium burning stars as a function of mean HI 
column density in square bins of width 1.7 arcmin 1.7 kpc). 
Error bars include Poisson noise and systematic error in the back- 
ground. Downward pointing arrows indicate those regions with 
only one young star. The line corresponds to a best fit slope of 
1.6 it 0.2 {2(t). Increasing or decreasing the area of each bin by a 
factor of 4 increases or decreases the slope by ~ 0.1 — 0.2, respec- 
tively. The measured slope is similar to the canonical disk-averaged 
Kenn icutt-Schmidt relation, which has a slope of 1.4 (jKennicutH 
UMI). 



imuthally averaged rings i n disk galaxies ( Wong fc Blitd 
I2002t iBoissier et al.l 120031 : iKennicuttI 120071 ). Recent the- 
oretical work on star formation in galactic disks predicts 
the KS relation to steepen at low gas densities, due to 
a star formation threshold, de creasing molecular frac- 
tion or flaring of the gas disk (iSchave fc Dalla Vecchial 
[2001 [Robertson fc Kravtsovl[2008f ). Since we have ne- 
glected molecular g detailed comparison to these 
studies is unwarranted. We simply note that our data 
are consistent with a single slope down to gas densities 
of 2 X 10^" cm~^, lower than typically probed using 
Ha flux as a star formation tracer and lower than most 
estimates of the star formation threshold, which include 
molecular gas. Note also that incompleteness may cause 
some of the fainter MS-t-BHeB stars to be preferentially 
missed in the most crowded regions with the highest gas 
densities, possibly biasing our measured slope downward. 
Nevertheless, it is encouraging that we recover a reason- 
able slope using an independent method sensitive to a 
wider age range than Ha flux (100 Myr vs. 10 Myr). 

5. RADIAL PROFILES 

We derived the V-band diffuse light surface brightness 
profile of M81 using the IRAF ellipse task with elliptical 
annuli of constant center, position angle, and inclination 
and after masking saturated stars. In each elliptical an- 
nulus, the median pixel value was computed after two 
5(7 clipping iterations. The sky value was estimated as 
the average median pixel value of 13 3x3 arcmin boxes 
near the edges of the mosaic. Squares in Fig. [T3] show 
the sky-subtracted profile. The error bars include read 
noise, Poisson noise in the total signal, and the stan- 
dard deviation of the 13 median sky values. Open circles 
represent the NE quadrant and filled circles represent 
the NW quadrant (relative to M81's nucleus and foUow- 
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Fig. 13. — V-band surface brightness profile of M81 derived from 
diffuse light (squares). Error bars are la and include read noise 
and uncertainties in the total signal and sky. Filled and open 
circles represent, respectively, the NW and NE quadrants following 
the convention where north is up and east is to the left of MSl's 
nucleus. The solid line is an exponential fit with a scale-length 
h = 2.7 ± 0.1 arcmin. The profile shows evidence for a fiattcning 
at ~ 17 arcmin. 



ing the convention where north is up and cast is to the 
left). The error bars on the circles are omitted for clarity. 
The median extinction of all point sources, Ay = 0.25, 
was subtracted from the profiles. Inside ~ 2 arcmin, the 
bulge is saturated, so there are no data points there. 

The diffuse light surface brightness profile traces MSl's 
exponential disk out t o Rdv /"^ 17 arcmin, farther 
than previou s studies dKmti [19871 : IMollenhofl [200l 
iWillner et al.l[200^ . The exponential scale-length inside 
17 arcmin (shown by the solid line) is 2.7 ± 0.1 arcmin, 
in good agreement with 2.53 arcmin found by MoUenhoff 
(2004), who performed a bulge-disk decomposition inside 
Rdp ~ 9 arcmin and over all 4 quadrants. The bump in 
the NE profile at ~ 9 arcmin is due to a spiral arm. 
There is a clear flattening in the profiles at ~ 17 arcmin, 
but the la uncertainties make the precise location and 
amount of flattening somewhat ambiguous. 

In Fig. [Ml we show the background-subtracted surface 
density of stars in each of the CMD selection boxes de- 
scribed in SJ31 The lines are color-coded so that the total 
RGB box is red, metal-poor RGB box is orange, metal- 
rich RGB box is green, AGB box is magenta, MS-t-BHeB 
box is blue, and RSG box is cyan. The top profile is 
the total of all the boxes. Each point in the profiles is 
the mean Rdp of all stars within a bin. We chose the 
bin widths so they would have approximately the same 
number of total stars. The vertical error bars include 
Poisson noise and systematic uncertainty in the back- 
ground, while the horizontal error bars denote the full 
radial range of stars in each bin. To count stars and 
pixel area, we used only those pixels with the highest 
weights in the confidence maps (i.e., with contributions 
from all dithered images). We also applied a small crowd- 
ing correction factor to the profiles, typ ically amounting 
to ^ 2 0%, following the prescription in llrwin fc Trimbig 
(|1984f ). This correction factor breaks down in the most 
crowded regions where the stellar PSFs significantly over- 
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Fig. 14. — Background-subtracted radial density profiles for stars 
in the different color-magnitude selection boxes. Vertical error 
bars include Poisson noise and systematic uncertainty in the back- 
ground. Horizontal error bars span the full radial range of stars in 
each bin. Severe crowding causes peaks in the RGB, AGB, and to- 
tal profiles at R^p ~ 10 — 14 arcmin. The total profile has a steep 
slope for i?£jp ~ 10 — 19 arcmin and a shallower slope for larger 
radii, confirming the behavior of the diffuse light profile. The ratio 
of metal-poor to metal-rich RGB star counts increases with radius, 
suggesting a decreasing mean metallicity. 

lap, resulting in a turnover in the density profiles at 
Rdp ~ 10 — 14 arcmin. Therefore, we limit most of the 
resolved star analysis to radii beyond this turnover. 

The young star profiles (MS-t-BHeB and RSG) exhibit 
the most substructure, as expected from their 2-D spa- 
tial distribution. The broad bump at 14 — 24 arcmin 
is due to HoIX and the most northern spiral arm. Be- 
yond 12 arcmin, the vast majority of the young stars 
are located in HoIX, Arp's Loop, and the northern two 
spiral arms. The sharp drop in the young star pro- 
files at Rdp ~ 10 — 12 arcmin approximately coincides 
with the star formation threshold r adius in M81 mea- 
sured bv lMartin fc KennicuttI (|2001h at Rhii « 12.3 ar- 
cmin. They obtained this value from the drop in the az- 
imuthally averaged Ha surface brightness profile over the 
whole disk, but asymmetries cause systematic azimuthal 
varia tions in Rhii , of a few arcmin llMartin fc Kennicutti 
I2OOII ). Like the young stars, the azimuthally averaged ra- 
dial profile of HI flux shows a peak at about 8 arcmin, 
with a pronounced hole near the nucleus and a steep 
gradient beyond 8 arcmin, so the peak in the young star 
profiles at this radius may be real rather than an artifact 
of incompleteness. 

The total star count profile in Fig. [T4l confirms the be- 
havior of the diffuse light profile, with a steep inner slope 
between 10 and 19 arcmin and shallower outer slope be- 
yond 19 arcmin. The break radius lies somewhat farther 
out than inferred from the diffuse light profile, but this 
difference is not significant given the la uncertainties in 
the latter. The fact that this transition to a flatter ra- 
dial falloff is seen in both diffuse light and star count 
analyses means it is genuine and not due to background 
uncertainties or completeness variations. 

The RGB stars dominate the total number counts out- 
side Rdp ^ 12 arcmin, so the total RGB profile has a 
similar flattening as the total star profile. Interestingly, 



the metal-poor RGB component has a shallower slope 
than the metal-rich component. The appearance of such 
a trend could have a number of causes. First, the back- 
ground counts in the metal-poor RGB box would have 
to be underestimated by a factor of ~ 2.5, or ^ 2.0 
if the metal-rich RGB box background were simultane- 
ously overestimated by a factor of ~ 10. Since these 
factors are much larger than the differences between the 
NGC 4244 and ACS contamination estimates and their 
statistical uncertainties, we consider it unlikely that the 
observed population gradient is due to an error in the 
background subtraction. It could be conceivable that the 
completeness of the metal-rich RGB box has a stronger 
dependence on galactocentric radius than the metal-poor 
RGB box because the metal-rich stars are fainter in the 
V-band. However, correcting for such a difference would 
actually enhance the population gradient by making the 
metal-rich profile rise even more sharply to smaller radii. 
Another possibility is that the photometric errors could 
be scattering an increasing fraction of stars from the 
metal-poor RGB box to the metal-rich box as radius de- 
creases and the stellar crowding increases. This also is 
unlikely because the median errors change by < 0.05 mag 
for radii > 15 arcmin and this amount is much smaller 
than the widths of the GMD boxes. 

The veracity of the population gradient apparent in 
the RGB density profiles is supported by Fig. [151 which 
compares the CMDs for three different radial ranges (in- 
dicated in the top of each panel) , moving outwards in ra- 
dius and starting where the two RGB boxes have roughly 
the same surface density. The similarity between the 
faint envelopes of the CMDs testifies further to the sta- 
bility of the photometric errors and completeness rate 
over this radial range. Overplottcd on the CMDs are 
theoretical RGBs from Marigo ct al. (200^) for [M/H] = 
-1.3, -0.7, and -0.4 and an age of 10 Gyr. Most of the 
sources at Iq > 24 and 0.2 < {V - I)o < 1.0 are un- 
resolved background galaxies. Each CMD has roughly 
the same number of objects, so the area covered by each 
and, consequently, the number of background galaxies 
progressively increases with radius. 

There is a clear excess of sources in all three CMDs 
where we would expect to find an RGB population. The 
color distribution of these sources shifts to the blue as ra- 
dius increases. In the innermost CMD, these sources are 
distributed roughly evenly between the isochrones, but in 
the outermost CMD, the highest density of sources lies 
between the -1.3 and -0.7 isochrones, indicating a peak 
[M/H] ~ —1.0. A similar result holds, but shifted down- 
ward by ^ 0.2 dex if we use the Dartmouth isochrones 
(|Dotter et al.ll20C)7| ) with the same age. Splitting the dif- 
ference between the stellar models, we estimate the ex- 
tended component has a peak metallicity [M/H] ~ — 1.1 
in the range Rdp = 32 — 44 arcmin. If the RGB stars were 
as young as 2 Gyr, then their metallicities would shift 
upward by ~ 0.4 dex. Because of the age-metallicity de- 
generacy, we cannot completely rule out age effects, but 
we consider it unlikely that the observed trend is entirely 
due to an age spread, as no single metallicity can span 
the full color range of the total RGB selection box. 

As discussed in 331 the Subaru point source catalog 
contains some unresolved background galaxies and com- 
pleteness becomes an important issue for RGB metal- 
licities [M/H] > —0.7, raising the question of whether 
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(V-I), (V-I), (V-I), 

Fig. 15. — Color-magnitude diagrams of point sources in three different deprojected ra dial ranges as indic ated in the panels. Overplotted 
are theoretical red giant branches for an age of 10 Gyr and [M/H] = -1.3, -0.7, and -0.4 IjMarigo et a l. 2008). The diagrams have the same 
total number of point sources. Most of the sources at Iq > 24 and 0.2 ^ {V — I)o < 1.0 are unresolved background galaxies. The red giant 
branch color shifts to the blue as radius increases, suggesting a decreasing metallicity. The right panel indicates the extended component 
has a peak metallicity in the range —1.3 < [M/H] < —0.7. 
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Fig. 16. — Comparing the color-magnitude diagrams for the ACS 
and Subaru data in fields 4 and 5 in Fig. □ (fidp ~ 20-30 arcmin). 
Overplotted are theoretical red giant branches for an age of 10 Gyr 
and [M/H] = -1.3, -0.7, and -0.4 (Marigo ct al. 2008). The 50% 
completeness level in the Subaru CMD is estimated to lie at /q ^ 
25.0 for (y-/)o < 1.0 and Iq ~ 24.4 for 1.0 < (V-I)o < 2.0 (see [[3] 
for details). The two diagrams appear reasonably similar, bearing 
in mind that the Subaru data have a brighter limiting magnitude 
and more background galaxies, which lie mostly at Iq > 24 and 
0.2 < (y - /)() < 1.0. 



background contamination, incompleteness, or system- 
atic color-magnitude shifts are biasing our metallicity es- 
timate. We address this question by comparing in Fig. [12] 
the ACS and Subaru ClMDs for the two ACS background 
fields. The ClVIDs contain point sources from the same 
two regions on the sky (fields 4 and 5 in Fig. [T|) , but the 
sources are not cross-matched. Thus, the Subaru CIVID 
contains true stars and unresolved background galaxies 
which, as previously discussed, are expected to mainly in- 



habit a region blueward of the RGB. The ClVIDs appear 
reasonably similar, and there is no obvious sign in the 
ACS CIVID of any significant metal-rich RGB population 
at this distance {Rdp ~ 20 — 30 arcmin) that could have 
been missed by the Subaru catalog. Based on the above 
considerations, we assign an uncertainty of ±0.3 dex to 
our metallicity estimate for the extended component. 

Usi ng RGB star cou n ts in s everal WFPC2 fields across 
IvISl, iTikhonov et al.l (|2005[ ) presented evidence for a 
flattening in stellar surface density at R^p ~ 20 arcmin. 
Their study, however, relied on only two fields beyond 
this radius on opposite sides of the galaxy (field 1 in Fig.[l] 
and another coincident with field 3). From the color of 
the RGB, they derived [Fe/H] = —0.6 in four of their 
fields at Rdp < 14 arcmin and [Fe/H] = -0.77 in field 1 of 
Fig- El lying at Rdp — 24 arcmin. IVIouhcine et al. (2005) 
also analyzed field 1 and derived [Fe/H] = —0.9 ± 0.3. 
The differences between their metallicity estimates and 
ours are not surprising because their estimates come from 
regions with a dominant or non-negligible contribution 
from the bright optical disk, which is more metal-rich 
than the extended component. IVIoreover, the extended 
component itself may have a metallicity gradient, so our 
estimate strictly applies to the region Rdp = 32 — 44 
arcmin. 

6. DISCUSSION 

Taking a similar approach as llrwin "eFall (l2005l) . we 
construct a composite surface brightness profile by using 
the diffuse light and total star count profiles in the re- 
gions where they are each most reliable. That is, within 
the bright optical disk, where the effects of incomplete- 
ness are most severe, we use the diffuse light profile 
because it is insensitive to these effects. In the outer 
regions, where the sky background dominates the dif- 
fuse light, but where the completeness rate is the highest 
and varies the least, we use the star counts, which have 
a higher contrast over the background than the diffuse 



14 



Barker et al. 



60 
7S 



T3 



18 
20 

22 

24 

26 

28 



O diffuse light 
□ star counts 



10 



20 30 



40 



50 



Deprojected Radius (arcmin) 



Fig. 17. — Composite V-band surface brightness profile for M81 
created by merging the diffuse light profile (diamonds) with the 
total stellar surface density profi le (squares). Errors on the diffuse 
light are the same as in Fig. 1131 and errors on the star counts are 
from Fig. [14] 



light. The overlapping region allows us to bring the star 
counts onto the same absolute luminosity scale as the 
diffuse light. By merging the diffuse light and resolved 
star count profiles in this way, we can trace M81's sur- 
face brightness over a larger radial range and to fainter 
magnitudes than is possible with either profile alone. 

The resulting composite V-band light profile (Fig. [Tf)) 
shows even more clearly the presence of two regimes, a 
steep inner gradient for R^p ^18 arcmin and a shallower 
outer gradient at larger radii. As shown in f}5l the inner 
profile gradient follows an exponential with h — 2.7 ar- 
cmin. Fitting the total star counts beyond 20 arcmin 
yields h = 12.9 ± 0.9 arcmin and 7 = 2.0 ± 0.2. These 
fits maximize the extended component's contribution at 
small radii because they ignore the bright optical disk. 

In Fig. 1181 we show two possible 3-component decom- 
positions of the composite light profile. The point sym- 
bols and error bars are the same as in Fig. [2l To rep- 
resent the bulge region, we have added the profile of 
iMollenhofl (|200l . In doing so, we assume his model 
is a good description of the true V-band fight distribu- 
tion. His total light profile inside 2 arcmin is shown as 
circles and his Sersic bulge component (with n = 3.0 
and i?e — 0.75 arcmin) is the long-dashed line. The 
bright optical disk with scale-length h = 2.7 arcmin is the 
short-dashed line. Because of background uncertainties 
and the limited radial baseline, the total profile's shape 
parameters cannot be very well constrained in a simul- 
taneous 3-component fit. Instead, we begin by using the 
halos of the MW and M31 as benchmarks for compari- 
son under the assumption that the extended component 
is M81's halo. Later, we consider other possibilities such 
as a thick disk. 

Numerous studies over the years have found that the 

MW stellar halo has a peak [Fe/H] 1.6 (or [M/H] - 

— 1.3 for a typical [a/Fe] ~ 0.4), and a power-law volume 
density distribution with exponent 7 -|- 1 ~ 3.0 (see the 
recent reviews by Helmi [2008] and Geisler et al. [2007]). 
iNewberg fc Yannvl (|200 61 also found the distribution of 
halo F-turnofi^ stars in the SDSS followed a Hernquist 
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Fig. 18. — Exploring possible decompositions of MSl's V-band 
surface brightness profile. Diffuse light (diamonds) and star counts 
(squares) are supplemented with the total bulge plus disk model of 
MoUenhoff (2004) (circles) to represent the bulge region {Rdp ^ 2 
arcmin) w here our images become saturated. Errors are the same 
as in Fig. 1171 Broken lines show individual components: Sersic 
bulge model of MoUenhoff (2004) with n = 3.0 and = 0.75 
arcmin (long dashed), exponential inner disk with scale-length 2.7 
arcmin (dashed), and two spherically symmetric Hernquist halos 
with scale radii of 14.0 arcmin (dot-dashed) and 53.1 arcmin (dot- 
ted). These scale radii are rec ent estimat es for the h alos of the 
MW UNewberg fc Yannvl [200^) and M31 Hlbata et al.i,200 7). re- 
spectively, scaled to the distance of M81 (1 arcmin ^ 1 kpc). Two 
solid lines show the total profiles after adding the bulge, inner disk, 
and one Hernquist halo. 



profile with best-fit scale radius « 14 kpc ^ . 

Several recent studies have found evidence for an ex- 
tended hal o in M31 wit h met allicities in the range - 
0.7 to -1.5 (iReitzel fc"G uhathak urtall2002l: [Kalirai et all 
120061 |Chapman.eLalJ 12006; RichardsoneLaU |2009D 
their photometric study of the southern quadrant of M31, 
llbata et al.l (j2007l ) showed that this component could be 
fit by a Hernquist profile with = 53.1 ± 3.5 kpc or 
an exponential profile with scale-length h = 46.8 ± 5.6 
kpc out to 150 kpc. Fitting the minor axis profile in re- 
gions devoid of spatial substructures, they estimated a 
projected power-law exponent of 7 = 1.91 ± 0.12. 

Two Hernquist profiles with the same scale radii as the 
MW and M31 halos translated to M81's distance are rep- 
resented by the dot-dashed and dotted lines, respectively, 
in Fig. [THl To plot these profiles, we assumed M81's ex- 
tended component is spherically symmetric and multi- 
plied the deprojected radius in the appropriate formulae 
by a geometrical factor of ^^0.5 (1 -I- cos'^{i)) ~ 0.8, to 
approximately account for depro jecting the h alos with 
M81's inclination, i = 58 deg ((Regan fc Vog cl 1994). 
The two Hernquist profiles provide an adequate descrip- 
tion of the data, but are virtually indistinguishable over 
the observed radial range. Extending the survey out to 
R(ip ^ 50 arcmin would offer the best chance of rul- 
ing out one Hernquist model or the other. Integrat- 
ing the profiles out to 100 arcmin, M81's halo would 
contain 10 — 15% of its total V-band luminosity, or 

^ We recall that 41% of the half-mass radius an d 55% of 
the ( projected) half-light radius of a Hernquist profile IIHernQuistI 
11990) . 
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Fig. 19. — Background-subtracted azimuthal variation in stellar 
surface density for projected radii of 12 - 18 arcmin. Error bars 
have the same meaning as those in Fig. 1141 The solid and dotted 
lines mark the major and minor axes of M81, respectively. The 
metal-poor RGB stars have a flatter profile, as expected if some 
of them belong to MSl's halo or a more face-on or thicker disk 
component. 



Lv ~ 3 — 6 X 10^ Lq. Uncertainties in the diffuse light 
sky subtraction in the range Rj^p 14—18 arcmin trans- 
late to a systematic uncertainty of about ±50% in our 
luminosity estimate, since this was the region used to 
normalize the star counts to absolute V-band luminosity. 
Nevertheless, MSl's extended component appears to be 
more luminous than the MW and M31 halos, which have 
Lv ~ 10^ Lq (|Carnev et al.lll990l : Hbata et al.ll2007| ). 

As we have discussed, the total star counts beyond 
Rdp 18 arcmin are dominated by metal-poor RGB 
stars. If this metal-poor RGB component is MSl's halo, 
then it should display a more circular morphology pro- 
jected on the sky than the metal-rich RGB component. 
We test this possibility in Fig. 1191 where we examine 
the background-subtracted azimuthal variation in stellar 
surface density in a circular annulus around MSl's nu- 
cleus. To minimize the effects of crowding on the profiles 
and maximize the baseline in position angle, we include 
only stars with projected radii of 12 - IS arcmin. This 
range spans deprojected radii of about 15 - 30 arcmin 
in total, with the precise Rdp sampled depending on the 
position angle. The solid and dotted vertical lines mark 
MSl's major and minor axes, respectively. 

The profiles have a roughly sinusoidal shape, as ex- 
pected for a disk morphology with MSl's inclination and 
position angle. The profile of the metal-poor RGB com- 
ponent appears marginally flatter than the metal-rich 
component, suggesting it could contain a superposition 
of disk and halo populations, or a disk structure with a 
larger scale height or lower inclination than the bright 
optical disk. Applying a K-S test on the azimuthal dis- 
tributions of the metal-poor and metal-rich RGB stars 
(in the annulus) gives a probability of 0.1% that they 
are drawn from the same parent distribution. The fact 
that the maxima of the curves are offset from the ma- 
jor axis is probably caused by the northern spiral arms, 
which host strong star formation west of the major axis. 

One may wonder whether the most recent close en- 
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Fig. 20. — Same points as in Fig. 1181 The dash-dotted line showrs 
the best-fit r^/"* law with effective radius Re = 1.8 it 0.1 arcmin. 
The dashed line is an exponential disk with scale-length of 2.7 
arcmin and the solid line is the total model. T he r ^/^ law provides 
a poorer fit than the Hernquist profiles in Fig. 1181 



counters with MS2 and NGC 3077 caused any measur- 
able asymmetries in the density distribution of MSl's 
halo. These encounters supposedly occurred ~ 200 — 300 
Myr ago, which is much less than the halo crossing time 
of 1 Gyr, so there should not have been enough time for 
the halo to reach a new equilibrium. We have examined 
the stellar density profiles on both sides of MSI, find- 
ing similar trends as seen in Fig. [TH This fact suggests 
the encounters did not signficantly disturb the halo, but 
additional data are needed at other position angles and 
reaching larger radii to confirm this suggestion. A survey 
of MSl's southern region, particularly along the minor 
axis, would help further constrain the interaction history 
and any role it had in shaping the halo (and disk). 

Fig. [20] tests the possibility that the extended compo- 
nent is the faint extension of MSl's bulge. In this fig- 
ure, the dashed line is the exponential disk with h — 2.7 
arcmin, the dot-dashed line is our best-fitting r^/^ law 
with effective radius Re = 1.8 ± 0.1 arcmin and the solid 
line is the sum of both components. In this case, the 
bulge contributes ~ 60% of the total V-band luminos- 
ity. The data favor the Hernquist profiles over the r^/^ 
law because the latter has a slope that is too steep at 
large radii. The r^/^ law would also imply a significant 
metallicity gradient in MSl's bulge, as lKong et al.l ()2000f ) 
estimated [M/H] ^ 0.2 inside 2 arcmin. Fitting a Sersic 
law yields n — 4.7 ± 0.3 and i?e — 1-7 ± 0.1 arcmin, but 
does not significantly improve the fit quality. 

Rather than a halo or bulge, MSl's extended compo- 
nent may be some type of disk structure. Both metal- 
poor and metal-rich RGB stars may belong to a single 
perturbed disk component with a negative metallicity 
gradient and a change in intrinsic orientation and/or 
scale-length at Rdp ~ IS arcmin. Alternatively, the 
extended component may be more analag ous to the 
MW' s thick disk or M31's extended disk (Ibata et all 
120051) , with distinct dynamics and formation history from 
the dominant inner disk. The MW's thick disk has a 
scale-length roughly 40% larger than the thin disk, or 
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~ 3.6 kpc (lJuric et al.ll2008|) . it has a mean metallicity 
([Fe/H]) ~ -0.6 (iGilmor'eet al. 1995; Ro bin et al.|[l99a 
ISouijiran et~al1l2003t lAllende Pr icto ct a l.ll2006D. and it 
contains 15% of the total disk light (B user et al.1119991 : 
IChen et all 120011: ILarsen fc Humphreys 2003'). In com- 
parison, M81's extended component is less centrally con- 
centrated, has a slightly lower metallicity, but has a sim- 
ilar luminosity. We note that the MW's thick disk has 
been studied mainly within a few kpc of the solar cylin- 
der, so it is difficult to draw firm conclusions from this 
comparison. 

Regardless of the disk/spheroid nature of the ex- 
tended component, M81 seems to be a nearby example 
of the so-cal l ed aii ti-tru ncated galaxy clas s ident ified by 
lErwin et all (pOOSl ) and lPohlen fc Truiillol (l200l. Some 
of the anti-truncated envelopes found by lErwin et al.l 
([2005') had rounder isophotes than the inner disks, sug- 
gesting a transition to a spheroidal component fo r those 
particular galaxies, but iPohlen fc Truiillol (|2006D could 
not reliably measure the outer isophote shapes in their 
sample. It is also worth noting that the break in MSl's 
profile occurs at a larger radius than nearly all the anti- 
truncated galaxies in Pohlcn fc Trujillo (2006). 

The results presented here may also have important 
ramifications for studies of the star for mation history 
(SFH) in M81. For field 6 in Fig. [J iWilliams et all 
pOOai derived a SFH in which ~ 60% of stars formed 
prior to 8 Gyr a go. In light of our re sults, the stellar 
population in the Willi ams et all ()200 8^ field is actually 
a mixture of the inner disk and extended component, 
which at this distance on the major axis could contribute 
~ 25 — 35% to the total V-band luminosity. We sug- 
gest, therefore, that if the extended component is an old 
spheroid or thick disk, then it could be responsible for 
some of the star formation in their oldest age bin (> 8 
Gyr), meaning M81's outer disk is younger than inferred 
from the total SFH. More generally, these considerations 
highlight the complications that can arise in pencil-beam 
studies when multiple galactic structural components lie 
along the same line-of-sight. Studies of the outer disks of 
massive spirals may contain a non-negligible "contamina- 
tion" from their spheroids or thick disks and vice-versa. 

7. SUMMARY 

Using Suprime-Cam on the Subaru telescope, we have 
conducted the first wide field mapping of both young 
(< 200 Myr) and old (~ 1 - 10 Gyr) resolved stellar 
populations around M81, reaching projected distances of 
~ 30 kpc. Throughout the surveyed region, the surface 
density of young stars (< 100 Myr) traces the HI column 
density in a manner similar to the Kennicutt- Schmidt 
relation, but we are able to probe the relation over a 
much wider range of gas densities than is normally pos- 
sible with Ha imaging. We find no strong evidence for a 
population of RGB stars in Arp's Loop over and above 
the background M81 population. This is consistent with 
the idea that Arp's Loop was initially completely gaseous 
tidal debris that began forming stars around the time of 
the last significant gravitational interaction of the M81 
group ~ 200 - 300 Myr ago. 

Both diffuse light and resolved star counts show evi- 
dence for an extended structural component beyond the 
bright optical disk with a much flatter surface bright- 



ness profile, but the star counts allow us to probe this 
component to fainter levels than is possible with the dif- 
fuse light alone. This component begins to dominate at 
Rdp ^ 18 arcmin and fiy ^ 26 mag arcsec^^, and con- 
tinues out to the last measured point at Rdp 37 arcmin 
and /iy ~ 28 mag arcsec"^. This flattening is accompa- 
nied by a shift in RGB color toward the blue, which we 
interpret as a shift to lower metallicity. 

The extended component shares some similarities with 
the MW's halo and thick disk, but has some important 
differences, as well. If this extended component is ^ 10 
Gyr old, then it has a peak metallicity [M/H] ^ — 1.1±0.3 
at radii Rdp = 32 — 44 kpc assuming a distance of 3.6 
Mpc. This metallicity is slightly higher than the MW's 
halo, but lower than the MW's thick disk. In the radial 
and azimuthal range probed by our data, the surface den- 
sity profile of the extended component follows a power- 
law with an exponent of 7 ~ 2, similar to the MW's 
halo, but much shallower than the MW's thick disk. If 
it is separate from M81's bulge and thin disk, then it 
contains ~ 10 — 15% of M81's total V-band luminos- 
ity, which is similar to the MW's thick disk, but several 
times more luminous than the MW's halo. Caution must 
be exercised when making such comparisons, however, 
since these structures have been isolated with different 
selection techniques and probed at different galactocen- 
tric distances. 

Other interpretations for the extended component, 
such as a perturbed thin disk, remain viable. An inter- 
action induced disturbance occurring < 1 Gyr ago seems 
an unlikely cause for this component, as it appears in 
both the NE and NW quadrants relative to the galaxy's 
nucleus. If this component is the faint extension of the 
bulge, a possibility which our data disfavor relative to a 
halo or thick disk, then the bulge would contain over half 
of the total V-band luminosity and must have a signifi- 
cant metallicity gradient, from super-solar in the nucleus 
to about 10 times less than solar at a projected radius of 
- 30 kpc. 

Further progress in understanding this structure in 
M81 requires mapping the resolved RGB stars over a 
wider area. This includes extending the radial baseline 
to properly test different structural models, and prob- 
ing more position angles to check for asymmetries that 
may have occurred as a result of gravitational interac- 
tions with its galactic neighbors. Kinematic information 
for individual RGB stars in this component would also 
be important for constraining its nature, but this must 
await future generations of ground-based telescopes and 
spectrographs. 
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